Solution of the Radial Schrodinger Equation for the Potential 
Family V{r) = 4~- + Cr K using the Asymptotic Iteration Method 

M. Aygunf, 0. Bayrak|§ and I. Boztosun§ 
^Faculty of Arts and Sciences, Department of Physics, 

Ataturk University, Erzurum, Turkey 
^Faculty of Arts and Sciences, Department of Physics, 

Bozok University, Yozgat, Turkey and 
^Faculty of Arts and Sciences, Department of Physics, 
Erciyes University, Kayseri, Turkey 
(Dated: February 3, 2008) 

Abstract 

We present the exact and iterative solutions of the radial Schrodinger equation for a class of 
potential, V(r) = 4~y + Cr K , for various values of k from -2 to 2, for any n and I quantum states 
by applying the asymptotic iteration method. The global analysis of this potential family by using 
the asymptotic iteration method results in exact analytical solutions for the values of k = 0,-1 
and —2. Nevertheless, there are no analytical solutions for the cases k = 1 and 2. Therefore, 
the energy eigenvalues are obtained numerically. Our results are in excellent agreement with the 
previous works. 
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I. INTRODUCTION 



We search for the exact solution of a new class of potential as in the following form: 

V(r) = --- + Cr K (1) 



This potential includes various potentials according to the values of the potential parameters: 
a) By choosing C = 0, this potential turns into the Kratzer potential, which is analytically 
solvable by using various methods . This potential has been used extensively so far 



flflQ. 



in order to describe the molecular structure and interactions 

b) The Gol'dman-Krivchenkov potential is obtained by choosing A ^ 0, B = and k = 2 
and it is solved analytically by using various methods 2J, |4j. Furthermore, this potential 
becomes the spiked harmonic oscillator potential 4] when A = 0, B = 1 and k — 2. 

c) The Coulomb plus linear potential form is obtained by choosing_A = and k — 1 and 
is solved by methods such as the envelope and variational methods [4j. This potential has 
raised great interest in atomic and molecular physics (Ref. j^] and references therein). 

d) The Coulomb plus harmonic oscillator potential form is obtained by choosing A = 
and k — 2. This form is solved by using the moment method the shifted 1/N expansion 
method [7( and the envelope method 4\. In addition, it is applied to examine the Zeeman 
quadratic effect 8] and the magnetic field effect in the hydrogen atom 

These potentials have been extensively used to describe the bound and the continuum 
states of the interactions systems and a great number of papers have been published for 
the exact and numerical solutions of these potentials (see Refs. 
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1 1 31 ] and references therein). Thus, it would be interesting and important to 
solve the non-relativistic radial Schrodinger equation for a new solvable potential family 
for any n and I quantum states. Recently, an alternative method, called as the Asymptotic 
Iteration Method (AIM), has been developed by Ciftgi et al. 13j for solving the second-order 
homogeneous linear differential equations and it has been applied to so 
radial Schrodinger equation as well as the relativistic wave equations 
In this paper, we aim to show that AIM could give the energy eigenvalues for the potentials 
that have analytical solutions obtained by using different methods. Moreover, AIM could 
provide energy eigenvalues for potentials that have no analytical solutions for any n and I 
quantum states with various k values. 



ve the non-relativistic 
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In the next section, we introduce AIM and then in section 1111} the analytical solution of 
the Schrodinger equation is obtained by using AIM for the potentials with k = 0,-1 and 
—2 and a closed form for the energy eigenvalues and corresponding eigenf unctions for any 
n and I quantum numbers are given. Then, in section IIVI the energy eigenvalues for a new 
solvable potential with k = 1 and 2 cases by using AIM iteration procedure are obtained. 
Finally, in section |Vj we remark on these results. 



II. OVERVIEW OF THE ASYMPTOTIC ITERATION METHOD (AIM) 



AIM is proposed to solve the second-order differential equations of the form 13j, Il6l |. 



y" = X (x)y' + s (x)y 



(2) 



where A (a;) ^ 0. The variables, s (x) and A (a;), are sufficiently differentiable. The differ- 
ential equation ([2]) has a general solution 



y[x) = exp 



a(x )dx 



(3) 



if k > 0, for sufficiently large k, we obtain the a(x) values from 

Sk(x) Sfc-l(x) 



X k (x) A fc -i(x) 



a(x), k = 1, 2, 3, 



(4) 



where 



X k (x) = Aj b _ 1 (x) + s k -i(x) + A (x)A fc _i(x) 

s k{x) = + s (x)X k ^ 1 (x), k = 1,2,3,. 



(5) 



It should be noted that one can also start the recurrence relations from k = with the 
initial conditions A_i = 1 and s_j = For a given potential, the radial Schrodinger 
equation is converted to the form of equation (|2J). Then, s (x) and X (x) are determined and 
Sfe(x) and X k (x) parameters are calculated by the recurrence relations given by equation (111). 

The energy eigenvalues are obtained from the roots of the quantization condition, given 
by the termination condition of the method in equation The quantization condition of 
the method together with equation (jSj) can also be written as follows 



= A fc (x)s fe ^i(a;) - A fe ^i(x)s fc (x) = k = 1,2, 3, 



(6) 



The energy eigenvalues are obtained from this equation if the problem is exactly solvable. If 
not, for a specific n principal quantum number, we choose a suitable xq point, determined 
generally as the maximum value of the asymptotic wave function or the minimum value of 
the potential, and the approximate energy eigenvalues are obtained from the roots of this 
equation for sufficiently great values of k with iteration. 

The wave functions are determined by using the following wave function generator 

y n {x) = C 2 exp{- J ^(x') dx '^ ^ 

where k > n,n represents the radial quantum number and k shows the iteration number. For 
exactly solvable potentials, the radial quantum number n is equal to the iteration number 
k and the eigenfunctions are obtained directly from equation ([7]). For nontrivial potentials 
that have no exact solutions, k is always greater than n in these numerical solutions and the 
approximate energy eigenvalues are obtained from the roots of equation ([6]) for sufficiently 
great values of k by iteration. 



III. n = 0, -1 AND -2 CASES: ANALYTICAL SOLUTIONS 



Inserting the potential given by equation ([!]) into the Schrodinger equation gives 



d 2 R n i 2m 
dr 2 h 2 



E 



A + B_ CrK _ l{l + l)h 2 



R 



nl 







(8) 



r' r 2mr 2 

where n and I are radial and orbital angular momentum quantum numbers, A, B and C are 
strictly positive constants. By using the following ansatze: 
o 2mE o ~ ~ 2mA 



— e 



h 2 



e nl 



t nl 



C, A 



h 2 



~ 2mB ~ 2mC 
B = — — — , C 



Equation (IE]) becomes 



d 2 R, 



nl 



dr 2 



+ 



— e 



A B ~ 1(1 + 1) 
H Cr K - 



n»2 Ijft fyt 2 



h 2 ' 



Rnl — 



h 2 



(9) 



(10) 



The aim of this paper is to show how to obtain the analytical and numerical solutions of 
V(r) = A — — + Cr K potential with different k values, k—0, -1 and -2 values of this potential 

n n n n 

have analytical solutions and have been extensively studied so far 121, 121, H8|]. Therefore, 
we do not give the details of the calculations and only show the analytical solutions. By 
inserting k — 0,-1 and —2 into equation ffTUj) and by doing some simple algebra described 



in the previous section, the energy eigenvalues and regular eigenfunctions are obtained as 
follows: 



k Eigenvalues Eigenfunctions 



0, Era = C - ^ (n + ± + J (I + \) 2 + *2£) , R nl (r) = Ne~<T°r lFl (-n, 2A K=0 + 2; 2e?=°* 



-1, £ ni = -^fg^ (n + i + ,/(/ + i) 2 + %^)~ 2 , = iVe-^r- 1 - ^(-n, 2A*=" 1 + 2; 2e 



2, = -^r [n+ ± + J(l + ±) 2 + p(A + C)) , ^(r) = Ne'^T - lFl (-n, 2A*=" 2 + 2; 2e K nl 



with 



,K=0 




fc ni 






1 _ 


F K= ~ 
t nl 






2 _ 


F K= ~ 
t nl 





= iMr r a— = -i + yJ(i + fr + A (ii) 

_ B Ak=-2 _ 1 



A .=-2 = _i + n l+ u2 + A + c 



2(ra+A+l)' ~'~ 2 

lent agreement with the previous results obtained by using different 



These results are in excei 

methods (see [l, 3, [l8] and references therein) 



IV. k = 1 AND 2 CASES: ITERATIVE SOLUTIONS 

For 1 and 2, the exact analytical solutions can not be found and in this 

section, we present how to find the energy eigenvalues by applying the asymptotic iteration 
method. If we consider V(r) = 4r — — + Cr K potential in the three-dimensional radial 
Schrodinger equation, we obtain equation (jHJ). The straightforward application of AIM to 
solve this equation gives us the energy eigenvalues, however, we have observed that the 
energy eigenvalues oscillate and do not converge within a reasonable number of iteration. 
The sequence appears to converge when the number of iterations k ~ 30, but then it begins 
to oscillate as the iteration number k increases. This result violates the principle behind 
the AIM; as the number of iteration increases, the method should converge and should not 
oscillate. In order to overcome this problem and obtain a rapid convergence, we make a 
change of variables as r = rop, where = 2~g> then we obtain 



d 2 R 



dp 2 

with the following ansatz 



p p 2 



Rnl = (12) 



H 2 E 2 2mCV*+ 2 . 1 /~ 1~~ ~ ~ 2mA . 

7 2 = ^^, ^-o + A/a+o) 2 + A A=— (13) 



2m5 2 ' ' h 2 ' 2 V v T ' h 2 



we can transform equation (fl2l) to another Schrodinger equation form by changing the 
variable to p = u 2 and then by inserting R(u) = u^(f){u) into the transformed equation. 
Therefore, we obtain the Schrodinger equation as follows 



d 2 (f>(u) 
du 2 



+ 



Aeu + 4 - 47% 



2„,2k+2 



A(A 



u- 



where A = 2l' + ± 



1 Case 







(14) 



For k — 1, we obtain the following equation by using equation fll4|) : 



u 



du 2 



+ 



A 2 A A 2 4 A(A+1) 

4ew 2 + 4 - 47 V i— ^ — - 



(j){u) = 



(15) 



If we take the wave function in the following form 



(u) = u A+1 exp(-^-)f(u) 



(16) 



where j3 is an arbitrarily introduced constant to improve the convergence speed of the method 
We also introduce 7 in the asymptotic wave function since the parameter 7 appears 
explicitly in the potential. It is known that when changing the value of 7 in the potential, it 
affects the eigenvalues and the shape of the wave function. Therefore, by introducing 7 into 
equation fflBT) . we control the change of the potential parameters and, as a result, obtain a 
better convergence. If we insert this wave function into equation (fl5|) . then we obtain the 
second-order homogeneous linear differential equation as follows; 



d 2 f{u) 
du 2 



2(2/? 7 u 3 - 



A + 1. 



u 



df(u) 



du 



+ [(4/37A + 10/37 - 4s)u 2 - 4/3 2 7 V - 4 + 4 7 V] f(u) 

(17) 



which is now amenable to an AIM solution. By comparing this equation with equation 
we can write the sq(u) and \q(u) values as below 



s (u) 
X (u) 



(4/5 7 A + 10/?7 - 4e)u 2 - 4/3 2 7 V - 4 + 4 7 V 
.3 A + 1, 



2(2/3 7 M d 



u 



-) 



(18) 
(19) 



In order to obtain the energy eigenvalues from equation flTjl by using equation ([5]), we 
obtain the Sfc(r) and A^(r) in terms of s (r) and A (r). Then, by using the quantization 



condition of the method given by equation (T4J), we obtain the energy eigenvalues. Therefore, 
we have to choose a suitable uq point to solve the equation S n (uo, e) = iteratively in order 
to find e values. In this study, we obtain the Uq from the maximum point of the asymptotic 

r \ 1/4 

wave function, which is the same as the root of Ao(u) = 0, thus u = i^p) • This 
straightforward application of AIM gives us the energy eigenvalues as the sequence appears 
to converge when the number of iterations (k) are k ~ 30 as shown in Table [IJ 

In this table, we also present the convergence rate of AIM calculations. The energy 
eigenvalues appear as the number of iterations are k ~ 30. However, the speed of the 
convergence depends on the arbitrarily introduced constant, f3. We have investigated the 
optimum values of f3 that give the best convergence and have kept the one that appears to 
yield the best convergence rate. Therefore, for case k=1, we have performed calculations for 
the different values of (3. It is seen from Table |T] that the best convergence is obtained when 
the constant (3 values are (3=0A, 0.5, and 0.6. For other values of /3, the convergence needs 
more iteration. 

Having determined the value of f3 in an empirical way, for case k—1, the energy eigenvalues 
by using AIM are shown for different values of n and I for convergence constant (3 = 0.5, 
A — 1 and 7 = 1 in Table [Til 



B. k = 2 Case 



For k = 2, we obtain the following equation by using equation ( Ti~4"l) : 
d 2 <b(u) 



du 2 



+ 



' 2 a a 2 6 A(A+1) 

Aeu 1 + 4 - 47 V - 



if 



(t*) = 



If we take the wave function in the following form 



u) = u A+1 exp(— )f(u) 



~f(3u 4 



(20) 



(21) 



When we insert this wave function into equation ( l20l . we obtain the second-order homoge- 
neous linear differential equation as follows 



d 2 f(u) 
du 2 



(2/? 7 u 3 - 



A + l, 



u 



df(u) 



du 



+ [(4/3 7 A + 10/37 - 4e)u 2 - 4/5 2 7 2 M 6 - 4 + 4 7 V] f(u) 

(22) 



which is now amenable to an AIM solution. By comparing this equation with equation 
we obtain sq(u) and \q(u) values as below 

s = 4/^7M 2 A - 4/3 2 7 V + 10/37M 2 - Aeu 2 - 4 + 4 7 V 



(23) 



A = 2(2/37^ - — ) (24) 
u 



Similar to case k — 1, in order to obtain the energy eigenvalues from equation (1221) by 
using equation ([5]), we obtain the Sfc(r) and Afc(r) in terms of s${r) and Ao(r). Then, by 
using the quantization condition of the method given by equation (TJJ, we obtain the energy 
eigenvalues. For this case, we choose u = and we obtain the energy eigenvalues after 

the k = 30 iterations. The AIM results are presented in Table II III for different values of 
n and I quantum numbers for different values of 7. In Table |IV] and |V] we compare our 
results with previous works conducted for special forms of our potential. Our results are 
in excellent agreement with moment method 6j, which is the closest one to exact solution. 
They are better than perturbation results. 

V. CONCLUSION 



In this paper, we have presented the exact and iterative solutions of the radial Schrodinger 
equation for a class of potential, V(r) = ^ — — + CV K , for various values of k from -2 to 2 for 
any n and I quantum states by applying the asymptotic iteration method. According to the 
value of k, this potential family includes Kratzer, Modified Kratzer, Goldman-Krivchenkov 
or spiked harmonic oscillator, Coulomb plus linear and Coulomb plus harmonic oscillator 
potentials. 

These potentials have been extensively used to describe the bound and the continuum 
states of the interactions systems and a great number of papers have been published for the 
exact and numerical solutions of these potentials (see Refs. 
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and references therein). We have examined these potentials in this paper and have 
attempted to obtain exact or numerical solutions. 

The global analysis of this potential family by using the asymptotic iteration method has 
resulted in exact analytical solutions for the values of k — 0, — 1 and —2. In addition, closed- 
forms for the energy eigenvalues as well as the corresponding eigenfunctions are obtained, 
but for k = 1 and 2, there are no analytical solutions. Therefore, the energy eigenvalues are 
obtained numerically. 

The advantage of the asymptotic iteration method is that it gives the eigenvalues directly 
by transforming the second-order differential equation into a form of y" =Ao(r)y / + so(r)y. 



S 



The wave functions are easily constructed by iterating the values of sq(t) and Xo(r). The 
asymptotic iteration method results in exact analytical solutions if there is and provides the 
closed-forms for the energy eigenvalues as well as the corresponding eigenfunctions. Where 
;here is no such a solution, the energy eigenvalues are obtained by using an iterative approach 
Jj], 16, 17, Jj], 19]- As it is presented, AIM puts no constraint on the potential parameter 



values involved and it is easy to implement. The results are sufficiently accurate for practical 
purposes. It is worth extending this method to examine other interacting systems. 
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2 36071 1 58 


2 36071 387 


2 36071 387 


2 36080271 
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2.36071418 


2.36076466 



TABLE I: For k = 1 case, the energy eigenvalues (e) for the n = and I = states by means of 
several (3 values. We take A = 1 and 7 = 1. 
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n = 0,1 = 


n = 1,Z = 


n = 1,1 = 1 


n = 2,1 = 


n = 2,1 = 1 
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6.073324 


6.704883 
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TABLE II: For k = 1 case, the energy eigenvalues (e) for the several quantum states. We take 
A = 1, 7 = 1 and convergence constant /3 = 0.5. 
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0.9621813797 


10.8379872414 


112.1314195487 


2 





0.9982832867 


11.5427585197 


120.1879082158 
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1.1684333950 


13.0102255287 


134.1850120876 




2 


1.372953445 


14.8691668817 


152.2273430525 



TABLE III: For k = 2 case, the energy eigenvalues (e) for several quantum numbers. We take 
A = 1 and convergence constant = 1. 

n I B C Shifted 1/JV E [7] Moment £[6] AIM £ [14] Present Work £ 

1 1 0.60025 0.59377 0.59365 0.59377 



TABLE IV: Comparison of our results with the shifted 1/JV expansion, moment and perturbative 
asymptotic iteration methods for the energy eigenvalues, where A = 0, B = 1, C = 1, k = 2, m=l 
and h = 1. 
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1000.0 
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59.37546904 


2000.0 


85.7348038 


85.73480386 


5000.0 


138.5571975 


138.55719764 



TABLE V: Comparison of our results with the ground state energy eigenvalues of the hydrogen 
atom in the V(r) = Cr 2 potential, where A = 0, B = 1, k = 2, m=l and h = 1. 
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